Estimating the population effectiveness of interventions against COVID-19 in France: A modelling study.

BACKGROUND: Non-pharmaceutical interventions (NPIs) and vaccines have been widely used to manage the COVID-19 pandemic. However, uncertainty persists regarding the effectiveness of these interventions due to data quality issues, methodological challenges, and differing contextual factors. Accurate estimation of their effects is crucial for future epidemic preparedness. METHODS: To address this, we developed a population-based mechanistic model that includes the impact of NPIs and vaccines on SARS-CoV-2 transmission and hospitalization rates. Our statistical approach estimated all parameters in one step, accurately propagating uncertainty. We fitted the model to comprehensive epidemiological data in France from March 2020 to October 2021. With the same model, we simulated scenarios of vaccine rollout. RESULTS: The first lockdown was the most effective, reducing transmission by 84 % (95 % confidence interval (CI) 83-85). Subsequent lockdowns had diminished effectiveness (reduction of 74 % (69-77) and 11 % (9-18), respectively). A 6 pm curfew was more effective than one at 8 pm (68 % (66-69) vs. 48 % (45-49) reduction), while school closures reduced transmission by 15 % (12-18). In a scenario without vaccines before November 2021, we predicted 159,000 or 168 % (95 % prediction interval (PI) 70-315) more deaths and 1,488,000 or 300 % (133-492) more hospitalizations. If a vaccine had been available after 100 days, over 71,000 deaths (16,507-204,249) and 384,000 (88,579-1,020,386) hospitalizations could have been averted. CONCLUSION: Our results highlight the substantial impact of NPIs, including lockdowns and curfews, in controlling the COVID-19 pandemic. We also demonstrate the value of the 100 days objective of the Coalition for Epidemic Preparedness Innovations (CEPI) initiative for vaccine availability.

Prediction of long-term humoral response induced by the two-dose heterologous Ad26.ZEBOV, MVA-BN-Filo vaccine against Ebola.

The persistence of the long-term immune response induced by the heterologous Ad26.ZEBOV, MVA-BN-Filo two-dose vaccination regimen against Ebola has been investigated in several clinical trials. Longitudinal data on IgG-binding antibody concentrations were analyzed from 487 participants enrolled in six Phase I and Phase II clinical trials conducted by the EBOVAC1 and EBOVAC2 consortia. A model based on ordinary differential equations describing the dynamics of antibodies and short- and long-lived antibody-secreting cells (ASCs) was used to model the humoral response from 7 days after the second vaccination to a follow-up period of 2 years. Using a population-based approach, we first assessed the robustness of the model, which was originally estimated based on Phase I data, against all data. Then we assessed the longevity of the humoral response and identified factors that influence these dynamics. We estimated a half-life of the long-lived ASC of at least 15 years and found an influence of geographic region, sex, and age on the humoral response dynamics, with longer antibody persistence in Europeans and women and higher production of antibodies in younger participants.

Modeling the kinetics of the neutralizing antibody response against SARS-CoV-2 variants after several administrations of Bnt162b2.

Because SARS-CoV-2 constantly mutates to escape from the immune response, there is a reduction of neutralizing capacity of antibodies initially targeting the historical strain against emerging Variants of Concern (VoC)s. That is why the measure of the protection conferred by vaccination cannot solely rely on the antibody levels, but also requires to measure their neutralization capacity. Here we used a mathematical model to follow the humoral response in 26 individuals that received up to three vaccination doses of Bnt162b2 vaccine, and for whom both anti-S IgG and neutralization capacity was measured longitudinally against all main VoCs. Our model could identify two independent mechanisms that led to a marked increase in measured humoral response over the successive vaccination doses. In addition to the already known increase in IgG levels after each dose, we identified that the neutralization capacity was significantly increased after the third vaccine administration against all VoCs, despite large inter-individual variability. Consequently, the model projects that the mean duration of detectable neutralizing capacity against non-Omicron VoC is between 348 days (Beta variant, 95% Prediction Intervals PI [307; 389]) and 587 days (Alpha variant, 95% PI [537; 636]). Despite the low neutralization levels after three doses, the mean duration of detectable neutralizing capacity against Omicron variants varies between 173 days (BA.5 variant, 95% PI [142; 200]) and 256 days (BA.1 variant, 95% PI [227; 286]). Our model shows the benefit of incorporating the neutralization capacity in the follow-up of patients to better inform on their level of protection against the different SARS-CoV-2 variants. Trial registration: This clinical trial is registered with ClinicalTrials.gov, Trial IDs NCT04750720 and NCT05315583.

Impact of variants of concern on SARS-CoV-2 viral dynamics in non-human primates.

The impact of variants of concern (VoC) on SARS-CoV-2 viral dynamics remains poorly understood and essentially relies on observational studies subject to various sorts of biases. In contrast, experimental models of infection constitute a powerful model to perform controlled comparisons of the viral dynamics observed with VoC and better quantify how VoC escape from the immune response. Here we used molecular and infectious viral load of 78 cynomolgus macaques to characterize in detail the effects of VoC on viral dynamics. We first developed a mathematical model that recapitulate the observed dynamics, and we found that the best model describing the data assumed a rapid antigen-dependent stimulation of the immune response leading to a rapid reduction of viral infectivity. When compared with the historical variant, all VoC except beta were associated with an escape from this immune response, and this effect was particularly sensitive for delta and omicron variant (p<10-6 for both). Interestingly, delta variant was associated with a 1.8-fold increased viral production rate (p = 0.046), while conversely omicron variant was associated with a 14-fold reduction in viral production rate (p<10-6). During a natural infection, our models predict that delta variant is associated with a higher peak viral RNA than omicron variant (7.6 log10 copies/mL 95% CI 6.8-8 for delta; 5.6 log10 copies/mL 95% CI 4.8-6.3 for omicron) while having similar peak infectious titers (3.7 log10 PFU/mL 95% CI 2.4-4.6 for delta; 2.8 log10 PFU/mL 95% CI 1.9-3.8 for omicron). These results provide a detailed picture of the effects of VoC on total and infectious viral load and may help understand some differences observed in the patterns of viral transmission of these viruses.

Leveraging Contact Network Information in Clustered Randomized Studies of Contagion Processes.

In a randomized study, leveraging covariates related to the outcome (e.g. disease status) may produce less variable estimates of the effect of exposure. For contagion processes operating on a contact network, transmission can only occur through ties that connect affected and unaffected individuals; the outcome of such a process is known to depend intimately on the structure of the network. In this paper, we investigate the use of contact network features as efficiency covariates in exposure effect estimation. Using augmented generalized estimating equations (GEE), we estimate how gains in efficiency depend on the network structure and spread of the contagious agent or behavior. We apply this approach to simulated randomized trials using a stochastic compartmental contagion model on a collection of model-based contact networks and compare the bias, power, and variance of the estimated exposure effects using an assortment of network covariate adjustment strategies. We also demonstrate the use of network-augmented GEEs on a clustered randomized trial evaluating the effects of wastewater monitoring on COVID-19 cases in residential buildings at the the University of California San Diego.

Impact of non-pharmaceutical interventions, weather, vaccination, and variants on COVID-19 transmission across departments in France.

BACKGROUND: Multiple factors shape the temporal dynamics of the COVID-19 pandemic. Quantifying their relative contributions is key to guide future control strategies. Our objective was to disentangle the individual effects of non-pharmaceutical interventions (NPIs), weather, vaccination, and variants of concern (VOC) on local SARS-CoV-2 transmission. METHODS: We developed a log-linear model for the weekly reproduction number (R) of hospital admissions in 92 French metropolitan departments. We leveraged (i) the homogeneity in data collection and NPI definitions across departments, (ii) the spatial heterogeneity in the timing of NPIs, and (iii) an extensive observation period (14 months) covering different weather conditions, VOC proportions, and vaccine coverage levels. FINDINGS: Three lockdowns reduced R by 72.7% (95% CI 71.3-74.1), 70.4% (69.2-71.6) and 60.7% (56.4-64.5), respectively. Curfews implemented at 6/7 pm and 8/9 pm reduced R by 34.3% (27.9-40.2) and 18.9% (12.04-25.3), respectively. School closures reduced R by only 4.9% (2.0-7.8). We estimated that vaccination of the entire population would have reduced R by 71.7% (56.4-81.6), whereas the emergence of VOC (mainly Alpha during the study period) increased transmission by 44.6% (36.1-53.6) compared with the historical variant. Winter weather conditions (lower temperature and absolute humidity) increased R by 42.2% (37.3-47.3) compared to summer weather conditions. Additionally, we explored counterfactual scenarios (absence of VOC or vaccination) to assess their impact on hospital admissions. INTERPRETATION: Our study demonstrates the strong effectiveness of NPIs and vaccination and quantifies the role of weather while adjusting for other confounders. It highlights the importance of retrospective evaluation of interventions to inform future decision-making.

Neutrophil Activation and Immune Thrombosis Profiles Persist in Convalescent COVID-19.

PURPOSE: Following a severe COVID-19 infection, a proportion of individuals develop prolonged symptoms. We investigated the immunological dysfunction that underlies the persistence of symptoms months after the resolution of acute COVID-19. METHODS: We analyzed cytokines, cell phenotypes, SARS-CoV-2 spike-specific and neutralizing antibodies, and whole blood gene expression profiles in convalescent severe COVID-19 patients 1, 3, and 6 months following hospital discharge. RESULTS: We observed persistent abnormalities until month 6 marked by (i) high serum levels of monocyte/macrophage and endothelial activation markers, chemotaxis, and hematopoietic cytokines; (ii) a high frequency of central memory CD4+ and effector CD8+ T cells; (iii) a decrease in anti-SARS-CoV-2 spike and neutralizing antibodies; and (iv) an upregulation of genes related to platelet, neutrophil activation, erythrocytes, myeloid cell differentiation, and RUNX1 signaling. We identified a "core gene signature" associated with a history of thrombotic events, with upregulation of a set of genes involved in neutrophil activation, platelet, hematopoiesis, and blood coagulation. CONCLUSION: The lack of restoration of gene expression to a normal profile after up to 6 months of follow-up, even in asymptomatic patients who experienced severe COVID-19, signals the need to carefully extend their clinical follow-up and propose preventive measures.

Using a population-based Kalman estimator to model the COVID-19 epidemic in France: estimating associations between disease transmission and non-pharmaceutical interventions.

In response to the COVID-19 pandemic caused by SARS-CoV-2, governments have adopted a wide range of non-pharmaceutical interventions (NPI). These include stringent measures such as strict lockdowns, closing schools, bars and restaurants, curfews, and barrier gestures such as mask-wearing and social distancing. Deciphering the effectiveness of each NPI is critical to responding to future waves and outbreaks. To this end, we first develop a dynamic model of the French COVID-19 epidemics over a one-year period. We rely on a global extended Susceptible-Infectious-Recovered (SIR) mechanistic model of infection that includes a dynamic transmission rate over time. Multilevel data across French regions are integrated using random effects on the parameters of the mechanistic model, boosting statistical power by multiplying integrated observation series. We estimate the parameters using a new population-based statistical approach based on a Kalman filter, used for the first time in analysing real-world data. We then fit the estimated time-varying transmission rate using a regression model that depends on the NPIs while accounting for vaccination coverage, the occurrence of variants of concern (VoC), and seasonal weather conditions. We show that all NPIs considered have an independent significant association with transmission rates. In addition, we show a strong association between weather conditions that reduces transmission in summer, and we also estimate increased transmissibility of VoC.

Modelling the response to vaccine in non-human primates to define SARS-CoV-2 mechanistic correlates of protection.

The definition of correlates of protection is critical for the development of next-generation SARS-CoV-2 vaccine platforms. Here, we propose a model-based approach for identifying mechanistic correlates of protection based on mathematical modelling of viral dynamics and data mining of immunological markers. The application to three different studies in non-human primates evaluating SARS-CoV-2 vaccines based on CD40-targeting, two-component spike nanoparticle and mRNA 1273 identifies and quantifies two main mechanisms that are a decrease of rate of cell infection and an increase in clearance of infected cells. Inhibition of RBD binding to ACE2 appears to be a robust mechanistic correlate of protection across the three vaccine platforms although not capturing the whole biological vaccine effect. The model shows that RBD/ACE2 binding inhibition represents a strong mechanism of protection which required significant reduction in blocking potency to effectively compromise the control of viral replication.

Design, immunogenicity, and efficacy of a pan-sarbecovirus dendritic-cell targeting vaccine.

BACKGROUND: There is an urgent need of a new generation of vaccine that are able to enhance protection against SARS-CoV-2 and related variants of concern (VOC) and emerging coronaviruses. METHODS: We identified conserved T- and B-cell epitopes from Spike (S) and Nucleocapsid (N) highly homologous to 38 sarbecoviruses, including SARS-CoV-2 VOCs, to design a protein subunit vaccine targeting antigens to Dendritic Cells (DC) via CD40 surface receptor (CD40.CoV2). FINDINGS: CD40.CoV2 immunization elicited high levels of cross-neutralizing antibodies against SARS-CoV-2, VOCs, and SARS-CoV-1 in K18-hACE2 transgenic mice, associated with viral control and survival after SARS-CoV-2 challenge. A direct comparison of CD40.CoV2 with the mRNA BNT162b2 vaccine showed that the two vaccines were equally immunogenic in mice. We demonstrated the potency of CD40.CoV2 to recall in vitro human multi-epitope, functional, and cytotoxic SARS-CoV-2 S- and N-specific T-cell responses that are unaffected by VOC mutations and cross-reactive with SARS-CoV-1 and, to a lesser extent, MERS epitopes. INTERPRETATION: We report the immunogenicity and antiviral efficacy of the CD40.CoV2 vaccine in a preclinical model providing a framework for a pan-sarbecovirus vaccine. FUNDINGS: This work was supported by INSERM and the Investissements d'Avenir program, Vaccine Research Institute (VRI), managed by the ANR and the CARE project funded from the Innovative Medicines Initiative 2 Joint Undertaking (JU).

SAMBA: A novel method for fast automatic model building in nonlinear mixed-effects models.

The success of correctly identifying all the components of a nonlinear mixed-effects model is far from straightforward: it is a question of finding the best structural model, determining the type of relationship between covariates and individual parameters, detecting possible correlations between random effects, or also modeling residual errors. We present the Stochastic Approximation for Model Building Algorithm (SAMBA) procedure and show how this algorithm can be used to speed up this process of model building by identifying at each step how best to improve some of the model components. The principle of this algorithm basically consists in "learning something" about the "best model," even when a "poor model" is used to fit the data. A comparison study of the SAMBA procedure with Stepwise Covariate Modeling (SCM) and COnditional Sampling use for Stepwise Approach (COSSAC) show similar performances on several real data examples but with a much reduced computing time. This algorithm is now implemented in Monolix and in the R package Rsmlx.

Targeting SARS-CoV-2 receptor-binding domain to cells expressing CD40 improves protection to infection in convalescent macaques.

Achieving sufficient worldwide vaccination coverage against SARS-CoV-2 will require additional approaches to currently approved viral vector and mRNA vaccines. Subunit vaccines may have distinct advantages when immunizing vulnerable individuals, children and pregnant women. Here, we present a new generation of subunit vaccines targeting viral antigens to CD40-expressing antigen-presenting cells. We demonstrate that targeting the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein to CD40 (αCD40.RBD) induces significant levels of specific T and B cells, with long-term memory phenotypes, in a humanized mouse model. Additionally, we demonstrate that a single dose of the αCD40.RBD vaccine, injected without adjuvant, is sufficient to boost a rapid increase in neutralizing antibodies in convalescent non-human primates (NHPs) exposed six months previously to SARS-CoV-2. Vaccine-elicited antibodies cross-neutralize different SARS-CoV-2 variants, including D614G, B1.1.7 and to a lesser extent B1.351. Such vaccination significantly improves protection against a new high-dose virulent challenge versus that in non-vaccinated convalescent animals.

Between-group comparison of area under the curve in clinical trials with censored follow-up: Application to HIV therapeutic vaccines.

In clinical trials, longitudinal data are commonly analyzed and compared between groups using a single summary statistic such as area under the outcome versus time curve (AUC). However, incomplete data, arising from censoring due to a limit of detection or missing data, can bias these analyses. In this article, we present a statistical test based on splines-based mixed-model accounting for both the censoring and missingness mechanisms in the AUC estimation. Inferential properties of the proposed method were evaluated and compared to ad hoc approaches and to a non-parametric method through a simulation study based on two-armed trial where trajectories and the proportion of missing data were varied. Simulation results highlight that our approach has significant advantages over the other methods. A real working example from two HIV therapeutic vaccine trials is presented to illustrate the applicability of our approach.

A model for establishment, maintenance and reactivation of the immune response after vaccination against Ebola virus.

The 2014-2016 Ebola outbreak in West Africa has triggered accelerated development of several preventive vaccines against Ebola virus. Under the EBOVAC1 consortium, three phase I studies were carried out to assess safety and immunogenicity of a two-dose heterologous vaccination regimen developed by Janssen Vaccines and Prevention in collaboration with Bavarian Nordic. To describe the immune response induced by the two-dose heterologous vaccine regimen, we propose a mechanistic ODE based model, which takes into account the role of immunological memory. We perform identifiability and sensitivity analysis of the proposed model to establish which kind of biological data are ideally needed in order to accurately estimate parameters, and additionally, which of those are non-identifiable based on the available data. Antibody concentrations data from phase I studies have been used to calibrate the model and show its ability in reproducing the observed antibody dynamics. Together with other factors, the establishment of an effective and reactive immunological memory is of pivotal importance for several prophylactic vaccines. We show that introducing a memory compartment in our calibrated model allows to evaluate the magnitude of the immune response induced by a booster dose and its long-term persistence afterwards.

Machine Learning and Mechanistic Modeling for Prediction of Metastatic Relapse in Early-Stage Breast Cancer.

PURPOSE: For patients with early-stage breast cancer, predicting the risk of metastatic relapse is of crucial importance. Existing predictive models rely on agnostic survival analysis statistical tools (eg, Cox regression). Here we define and evaluate the predictive ability of a mechanistic model for time to distant metastatic relapse. METHODS: The data we used for our model consisted of 642 patients with 21 clinicopathologic variables. A mechanistic model was developed on the basis of two intrinsic mechanisms of metastatic progression: growth (parameter α) and dissemination (parameter µ). Population statistical distributions of the parameters were inferred using mixed-effects modeling. A random survival forest analysis was used to select a minimal set of five covariates with the best predictive power. These were further considered to individually predict the model parameters by using a backward selection approach. Predictive performances were compared with classic Cox regression and machine learning algorithms. RESULTS: The mechanistic model was able to accurately fit the data. Covariate analysis revealed statistically significant association of Ki67 expression with α (P = .001) and EGFR expression with µ (P = .009). The model achieved a c-index of 0.65 (95% CI, 0.60 to 0.71) in cross-validation and had predictive performance similar to that of random survival forest (95% CI, 0.66 to 0.69) and Cox regression (95% CI, 0.62 to 0.67) as well as machine learning classification algorithms. CONCLUSION: By providing informative estimates of the invisible metastatic burden at the time of diagnosis and forward simulations of metastatic growth, the proposed model could be used as a personalized prediction tool for routine management of patients with breast cancer.

Comparison of empirical and dynamic models for HIV viral load rebound after treatment interruption.

OBJECTIVE: To compare empirical and mechanistic modeling approaches for describing HIV-1 RNA viral load trajectories after antiretroviral treatment interruption and for identifying factors that predict features of viral rebound process. METHODS: We apply and compare two modeling approaches in analysis of data from 346 participants in six AIDS Clinical Trial Group studies. From each separate analysis, we identify predictors for viral set points and delay in rebound. Our empirical model postulates a parametric functional form whose parameters represent different features of the viral rebound process, such as rate of rise and viral load set point. The viral dynamics model augments standard HIV dynamics models-a class of mathematical models based on differential equations describing biological mechanisms-by including reactivation of latently infected cells and adaptive immune response. We use Monolix, which makes use of a Stochastic Approximation of the Expectation-Maximization algorithm, to fit non-linear mixed effects models incorporating observations that were below the assay limit of quantification. RESULTS: Among the 346 participants, the median age at treatment interruption was 42. Ninety-three percent of participants were male and sixty-five percent, white non-Hispanic. Both models provided a reasonable fit to the data and can accommodate atypical viral load trajectories. The median set points obtained from two approaches were similar: 4.44 log10 copies/mL from the empirical model and 4.59 log10 copies/mL from the viral dynamics model. Both models revealed that higher nadir CD4 cell counts and ART initiation during acute/recent phase were associated with lower viral set points and identified receiving a non-nucleoside reverse transcriptase inhibitor (NNRTI)-based pre-ATI regimen as a predictor for a delay in rebound. CONCLUSION: Although based on different sets of assumptions, both models lead to similar conclusions regarding features of viral rebound process.

Properties and pitfalls of weighting as an alternative to multilevel multiple imputation in cluster randomized trials with missing binary outcomes under covariate-dependent missingness.

The generalized estimating equation (GEE) approach can be used to analyze cluster randomized trial data to obtain population-averaged intervention effects. However, most cluster randomized trials have some missing outcome data and a GEE analysis of available data may be biased when outcome data are not missing completely at random. Although multilevel multiple imputation for GEE (MMI-GEE) has been widely used, alternative approaches such as weighted GEE are less common in practice. Using both simulations and a real data example, we evaluate the performance of inverse probability weighted GEE vs. MMI-GEE for binary outcomes. Simulated data are generated assuming a covariate-dependent missing data pattern across a range of missingness clustering (from none to high), where all covariates are measured at baseline and are fully observed (i.e. a type of missing-at-random mechanism). Two types of weights are estimated and used in the weighted GEE: (1) assuming no clustering of missingness (W-GEE) and (2) accounting for such clustering (CW-GEE). Results show that, even in settings with high missingness clustering, CW-GEE can lead to more bias and lower coverage than W-GEE, whereas W-GEE and MMI-GEE provide comparable results. W-GEE should be considered a viable strategy to account for missing outcomes in cluster randomized trials.

Dynamics of the Humoral Immune Response to a Prime-Boost Ebola Vaccine: Quantification and Sources of Variation.

The Ebola vaccine based on Ad26.ZEBOV/MVA-BN-Filo prime-boost regimens is being evaluated in multiple clinical trials. The long-term immune response to the vaccine is unknown, including factors associated with the response and variability around the response. We analyzed data from three phase 1 trials performed by the EBOVAC1 Consortium in four countries: the United Kingdom, Kenya, Tanzania, and Uganda. Participants were randomized into four groups based on the interval between prime and boost immunizations (28 or 56 days) and the sequence in which Ad26.ZEBOV and MVA-BN-Filo were administered. Consecutive enzyme-linked immunosorbent assay (ELISA) measurements of the IgG binding antibody concentrations against the Kikwit glycoprotein (GP) were available for 177 participants to assess the humoral immune response up to 1 year postprime. Using a mathematical model for the dynamics of the humoral response, from 7 days after the boost immunization up to 1 year after the prime immunization, we estimated the durability of the antibody response and the influence of different factors on the dynamics of the humoral response. Ordinary differential equations (ODEs) described the dynamics of antibody response and two populations of antibody-secreting cells (ASCs), short-lived (SL) and long-lived (LL). Parameters of the ODEs were estimated using a population approach. We estimated that half of the LL ASCs could persist for at least 5 years. The vaccine regimen significantly affected the SL ASCs and the antibody peak but not the long-term response. The LL ASC compartment dynamics differed significantly by geographic regions analyzed, with a higher long-term antibody persistence in European subjects. These differences could not be explained by the observed differences in cellular immune response.IMPORTANCE With no available licensed vaccines or therapies, the West African Ebola virus disease epidemic of 2014 to 2016 caused 11,310 deaths. Following this outbreak, the development of vaccines has been accelerated. Combining different vector-based vaccines as heterologous regimens could induce a durable immune response, assessed through antibody concentrations. Based on data from phase 1 trials in East Africa and Europe, the dynamics of the humoral immune response from 7 days after the boost immunization onwards were modeled to estimate the durability of the response and understand its variability. Antibody production is maintained by a population of long-lived cells. Estimation suggests that half of these cells can persist for at least 5 years in humans. Differences in prime-boost vaccine regimens affect only the short-term immune response. Geographical differences in long-lived cell dynamics were inferred, with higher long-term antibody concentrations induced in European participants.